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Abstract 



We compute the Renormalization Group functions of a Landau-Ginzburg- 
Wilson Hamiltonian with 0(n)xO(m) symmetry up to five- loop in Minimal 
Subtraction scheme. The line n~^{m,d), which limits the region of second- 
order phase transition, is reconstructed in the framework of the e = 4 — d 
expansion for generic values of m up to O(e^). For the physically interesting 
case of noncollinear but planar orderings (m = 2) we obtain n^{2, 3) = 6.1(6) 
by exploiting different resummation procedures. We substantiate this results 
re-analyzing six-loop fixed dimension series with pseudo-e expansion, obtain- 
ing ri+(2, 3) = 6.22(12). We also provide predictions for the critical exponents 
characterizing the second-order phase transition occurring for n > n'^. 

PACS Numbers: 05.70.Jk; 64.60.Pr; 75.10.Hk; ll.lO.Kk 
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I. INTRODUCTION 



In the last years much effort has been dedicated in the quahtative understanding and 
quantitative description of frustrated spin systems with noncolhnear order. Despite the 
intensive theoretical and experimental investigation, strongly debated issues remain open. 
In three dimensions it is not clear if systems as stacked triangular antiferromagnets (STAs) 
and Helimagnets should undergo a continuous phase transitions (named chiral) belonging 
to an unusual universality class and hence characterized by new values of critical exponents, 
or if the transition is first order (see Refs. [1-4] as review). 

Field-theoretical studies of such systems are based on the 0(n)xO(m) symmetric 
Landau- Ginzburg- Wilson (LGW) Hamiltonian [5,1] 




where (f)a o, < m) are m sets of n-component vectors. Whether m < n and Vq > 0, Eq. 
(1) describes m-dimensional spin orderings in isotropic n-spin space with the 0(n) — >0(n — 
m) pattern of symmetry breaking. Negative values of Vq correspond to simple ferromagnetic 
or antiferromagnetic ordering, and to magnets with sinusoidal spin structures [1]. The 
condition < vq < m/{m — 1)uq is required to have noncoUinearity and the boundedness 
of the free energy. For m — 2 the Hamiltonian (1) describes magnets with noncolhnear but 
planar orderings as frustrated XY {n — 2) and Heisenberg {n — 3) antiferromagnets, while 
for m > 3 systems with noncoplanar ground states [6,1]. 

The RG flow has always two flxed points (FP's), the Gaussian one (mq = fo = 0) and 
the 0{nm) one {vq = 0), while other two FP's (called chiral and antichiral) appear for some 
values of n. At fixed m, for the physically relevant case m > 1, the e expansion predicts 
four different regimes characterized as follows [7]: 

1) For n > n^{m, d), there are four FP's, and the chiral one is stable. 

2) For n~{m, d) < n < n'^{m, d), only the Gaussian and the Heisenberg 0(nm)-symmetric 
FP's are present, and none of them is stable. Thus the system is expected to undergo 
a first-order phase transition. 

3) For nH{rn,d) < n < n~{m,d), there are again four FP's. However for m < 7 the 
stable FP is in the region v < and the transition for those systems with > is 
again first order. 

4) For n < nnim, d), ^ there are again four FP's, and the Heisenberg 0(nm)-symmetric 
one is stable. 



^The precise value of n//(m, d) may be inferred from the results of Refs. [8,9], where it was shown 
that at a global 0{N) FP (here = nm) all the spin four operators (and the vq term is spin four) 
have the same scaling dimension. Thus the 0{nm) FP is stable for n < Nc/m = hh, where Nc is 
the marginal spin dimensionality of the cubic model that results in Nc ^ 2.9 in three dimensions. 
See the review [3] for an updated list of estimates and references about Nc- 
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For m > 3 all field theoretical (FT) methods (e = 4 — d [6,10], 1/n [10,11], fixed di- 
mensional d = ?> [12], e = d — 2 [13] expansions, and the Effective Average Action Method 
(EAAM) [14]) agree indicating the existence of a marginal number of components of the 
order parameter n^{m, 3) above which a second-order phase transition is driven by a stable 
chiral FP. 

For m = 2, again all the above mentioned FT methods indicate the existence of n^(2) ~ 6 
above which we expect a continuous phase transition [15,10,11,13,16-20,25-27]. Instead, for 
n < the situation is controversial. Several magnetic materials that are supposed to be 
described by the Hamiltonian (1) with n — 2,3 display scaling laws, but with experimental 
measured exponents that apparently depend on the material and sometimes also from the 
measure (see e.g. the recent experimental works [21] and also Refs. [3,4] for an exhaustive 
list of experiments). The same scenario is found in MC simulations: systems like STA dis- 
play scaling with almost well defined exponents [22] (see also [23] as the only simulation 
where apparently there is a first-order transition in STA), whereas systems with the same 
symmetry like rigid STA (STAR) and Stiefel Vn,2 (see e.g. Ref. [4] for the definition of 
these systems) clearly undergo first-order transition for n = 2, 3 [24]. Unfortunately a lot 
of theoretical efforts have not yet produced a unified picture. EAAM does not provide any 
stable FP favoring a weak first-order phase transition with nonuniversal pseudo-critical ex- 
ponents [25-27,4] that are in rough agreement with those coming from scaling relation found 
in experiments and Monte Carlo simulations. Whereas in six-loop fixed dimension approach 
another marginal number of components of the order parameter Nc2 < exists [20], so 
that for Nc2 < n < there is no stable FP, while for n < a stable (focus-like) chiral FP 
appears [18,19]. Since Nc2 — 5.7(3) [20], according to this scenario, XY and Heisenberg frus- 
trated antiferromagnets undergo a second-order phase transition, if they are in the domain 
of attraction of the stable FP. The found critical exponents [18,28] are in rough agreement 
with Monte Carlo and experimental results. Anyway, the continuous phase transition at this 
FP is very peculiar: scaling properties are governed over several decades of temperatures, by 
preasymptotic effective exponents, which can differ significantly from the asymptotic ones, 
explaining in this manner the apparently scattering experimental and numerical results [19]. 
Note that the existence of a stable FP docs not exclude the possibility that some systems 
may undergo a first-order transition. Indeed, they may lie outside the attraction domain 
of the stable FP. According to mean-field arguments the systems undergoing second-order 
phase transition are those characterized by Vq/uq < m/{m — 1). Since RG iterations may 
only narrow this region, systems that are outside this region surely undergo a first-order 
phase transition. 

The stable FP found in fixed dimension for n < Nc2 is unrelated to the small-e and large- 
n chiral one.^ This fact does not imply an inconsistency of all perturbative results, since e 



^This statement may seem at first sight far to be trivial. Anyway it is quite simple to understand 
for n = m = 2. In this case the e expansion has four FP's, but those with v* ^ Q are in the 
negative half plane (they are the counterpart of the tetragonal model according to the mapping of 
Rcf. [16]). Since in e expansion a theory with n couplings has at most 2" FP's, and the 0(2) x 0(2) 
has already four FP's no other FP may exist in e expansion. Thus for n = m = 2 it is impossible in 
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TABLE I. Values of n+(m, 3) for m = 2, 3, 4, 5 obtained with different FT methods. 





Method 


m=2 


m=3 


m=4 


m=5 


Ref. [32] 1994 


Local Potential Approximation 


~ 4.7 








Ref. [16] 1994 


d = 3 expansion: 0{g'^) 


3.91(1) 








Ref. [15] 1995 


e expansion: O(e^) 


3.39 








Ref. [25] 2000 


EAAM 


^ 4 








Ret. [27]2ooi 


EAAM 


~ 5 








Ref. [10]2ooi 


1/n expansion: 0(l/n^) 


5.3 


7.3 


9.2 


11.1 


Ref. [10]2ooi 


e expansion: O(e^) 


5.3(2) 


9.1(9) 


12(1) 




Ref. [11] 2002 


1/n expansion: 0(l/n^) 


~ 3.24 








Ref. [20] 2003 


d = 3 expansion: 0{g'^) 


6.4(4) 


11.1(6) 


14.7(8) 


18(1) 


This work 


€ expansion: O(e^) 


6.1(6) 


9.5(5) 


12.7(7) 


15.7(1.0) 


This work 


pseudo-e expansion: O(t^) 


6.22(12) 


9.9(3) 


13.2(6) 


16.3(1.3) 



and 1/n expansions describe adiabatic moving from 4 dimensions and n = oo respectively 
and they are probably inadequate to describe the essentially three-dimensional features of 
the chiral FP for physical systems. 

In recent years many efforts have been made to obtain a precise determination of 
n"'"(3, m). A complete list of results may be found in Table I. These scattered results 
clearly indicate that the extrapolation of low order 1/n and e expansion predictions up to 
physical values of n and e is a quite delicate matter. For this reason we extend the knowledge 
of the RG function up to five-loop in e expansion. Furthermore, to improve the goodness of 
fixed-dimensional = 3 estimates we re-analyze the six- loop perturbative series [18,12] with 
the pseudo-e expansion trick [29], since in many cases this method provided the most ac- 
curate results in the determination of the marginal number of order parameter components 
(see Ref. [30] for the cubic model and Ref. [31] for the weakly diluted m- vector model). We 
anticipate that also in this case the final estimate is very precise, as clear from Table I. 

The paper is organized as follows. In Section II we determine n'^(m, 3) for several n by 
means of five-loop e expansion. Pseudo-e expansion is considered in Section III. Section 
IV summarizes our results. In the Appendix A we report the estimates of the exponents 
for m = 2 and n > n"^(2, 3), whereas in the Appendix B details of series computation are 
presented, and the Appendix C reports analytic forms of some quantities given in the text. 



e expansion to find a signature of the three-dimensional chiral FP of Rcfs. [18,19]. For n = 3 and 
m = 2 the situation is different since there are two real FP's and two complex zeros of /3s. Anyway 
the latter have almost vanishing v coordinate, and they are far from the three-dimensional chiral 
FP, ruling out the idea that there is some connection between them. 
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II. FIVE-LOOP e-EXPANSION 



We extend the three-loop e expansion of Refs. [15,10] for the RG functions of the 
0(n)xO(m) symmetric theory to five-loop. The obtained series and some details of calcu- 
lations are reported in the Appendix B. Within these series, the e expansion of n^{m, 4 — e) 
may be calculated to O(e^). They are expanded as 



n^(m) = riQ (m) + nf{m) e + n^im) + nf{m) + n^{m) + O(e^), 
and the coefficients nf{rn) may be obtained by requiring 



and 



det 



d{u, v) 



(2) 



(3) 



(4) 



Note that the previous equation is not the only possible choice to find n^, one can e.g. 

impose the coincidence of the coordinates of the chiral and antichiral FP etc. However this 
was the most advantageous from the computational point of view. For generic values of m, 
the results obtained for the location of the FP's and for n^(m) are too cumbersome in order 
to be reported here, thus we report them only at fixed m = 2, 3. 



A. NoncoUinear but planar ordering: vii—2 

Let us first consider the physically relevant case m — 2. The numerical expression for 
n^(2,4 — e) simplifies to (their analytical expressions are in App. C) 

n+(2, 4 - e) = 21.7980 - 23.4310 e + 7.0882 - 0.0321 + 4.2650 + O(e^) , (5) 
n-(2, 4 - e) = 2.2020 - 0.5690 e + 0.9892 - 2.2786 + 6.5406 + O(e^) . 

The first three terms of the expansions are in agreement with previous works [15,10]. In order 
to give a numerical estimate of ri^{2, 3) such series should be evaluated at e = 1. A direct 
sum is obviously not effective since the series are clearly divergent. The series n~(2,4 — e) 
has a quite regular behavior in e, with the last terms approximately factorial growing and 
with alternating sign coefficients. Thus a Pade-Borel-Leroy (PEL) resummation (see App. 
B) should be effective, and in fact it provides n~(2,3) = 1.968(1). Despite to the stability 
of the result, this estimate cannot be right. Indeed we known, from the mapping onto the 
tetragonal model [16], that for n = m = 2 a couple of FP's exists in the e expansion, thus 
rr{2, 3) > 2 holds. A direct analysis of three-dimensional six-loop series shows that n~(2, 3) 
is a bit larger than 2 [33]. This was firstly noted in Ref. [16]. Anyway the precise value 
of n~(2,3) is not of interest for frustrated models and this disagreement will be not more 
commented below. We only quote it to show that PBL method apparently underestimates 
the final error bar. 
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TABLE II. Non-defective PBL results for n"'"(m = 2,3). The estimate from these data is 5.47(7). 



Loop 
Fade 


3 
[1/1] 


[1/2] 


4 


[2/1] 


[1/3] 


5 


[3/1] 


b=0 


3.385 


4.477 




5.423 


5.402 




5.455 


b=l 


3.514 


4.593 




5.423 


5.586 




5.455 


b=2 


3.583 


4.654 




5.423 


5.688 




5.455 


b=3 


3.626 


4.693 




5.423 


5.753 




5.455 



The evaluation of n+(2, 3) is much more difficult because of its irregular behavior. The 
results of PBL resummations are shown in Table II for several approximants. The fifth order 
ones indicate n+(2, 3) = 5.47(7) (it is the average and the variance of the approximants [1/3] 
and [3/1] with b = 0, 1). Anyway the previous analysis of n~(2, 3) suggests that PBL method 
underestimates error bars, thus to have a more reliable result, wc apply other summation 
methods to evaluate n+(2,3), that in practice allow one to extract better behaved series 
from Eq. (5). In what follows we present three different functions of n"'"(2,3) that behave 
better than Eq. (5). The choice of the functions is motivated or by physical reasons (as Eq. 
(7)) or by the experience in resumming divergent series (as Eq. (6)). Obviously we do not 
exclude the possibility of better choices. 

The first function (as usually done for critical exponents) is 

-—^ = 0.0459 + 0.0493e + 0.0381e^ + 0.0250e^ + 0.0056e^ + O(e^) , (6) 

n+(2,4 - e) v / > v y 

whose coefficients decrease rapidly. Setting e = 1, without resummation,^ one obtains 
n+(2,3) = 7.50 (at 3-loop), n+(2,3) = 6.31 (at 4-loop), and n+(2,3) = 6.10 (at 5-loop). 

Another method, firstly employed in Ref. [10], use the knowledge of n+(2, 2) to constrain 
the analysis at e = 2, under the assumption that n'^{m,d) is sufficiently smooth in d at 
fixed m. In Ref. [10] it is assumed that the two-dimensional LGW stable FP is equivalent 
to that of the NLcr model for all n > 2, m > 2, except n = 2, m = 2. Since the NLcr model 
is asymptotically free, the authors of Ref. [10] conclude that n^{2,2) = 2. The knowledge 
of n"'"(2,2) may be exploited in order to obtain some informations on n"'"(2,3), rewriting 
n+(2,4 - e) as 

n+(2, 4 - e) = 2 + (2 - e)(9.8990 - 6.7660e + 0.1611e2 + 0.0645e^ + 2.1648e^) + O(e^) 

= 2 + (2-e)a(m = 2,e). (7) 

Since the coefficients of a(2,e) do not decrease, we consider l/a(2,e) obtaining the more 
"convergent" expression 



^We do not use any resummation since the series is convergent up to the considered order. A 
good resummation should reproduce the result of the direct sum. Anyway note that PBL is not 
expected to work, since the series has not alternating signs. The same remark holds for Eq. (8). 
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n"'(2,3) 



^ - ^ method IV 

<$► ^ method III 

- ■ - ^ method II 
-©method I 

"1 r 

3 4 5 

Number of loops 

FIG. 1. Different estimates of n+(2,3) with varying the number of loops. Method I is PBL, 
method II is Eq. 6, method III is Eq. 8, and method IV is the same of III, but with n+(2,2) = 4. 
The shadow region is our final estimate n^(2, 3) = 6.1(6). 



— = 0.1010 + 0.0690e + 0.04566^ + 0.0294e3 - 0.0032e^ + O(e^), (8) 

ayZ, ej 

which, setting e = 1, gives n+(2,3) = 6.64 at three loop, ^+(2,3) = 6.08 at four loop, and 
n+(2,3) = 6.136 at five loop. 

Actually the exact value of n+(2,2) is a controversial matter, in fact it was pointed 
out that Z2-type topological defects may lead to a finite-temperature phase transition in 
the two-dimensional 0(3) x 0(2) model (see e.g. Refs. [34,35] and references therein, for 
different scenarios about the critical behavior of this model). Independently from the fact 
that the perturbative LGW approach is able to describe such topological excitations (that is 
still a very controversial point), one should expect as naive upper bound ^"''(2, 2) = 4, since 
for n > 4 the NLcx surely works. Using as constraint first ^^{2, 2) = 2 and then n"'"(2, 2) = 4 
one should kept in between the right value. Taking n"''(2,2) = 4 and considering again 
the series of the inverse of the resulting a(2,e), we have n^{2,3) = 7.57 at three loop, 
n~*"(2,3) = 6.91 at four loop, and n~^{2,3) = 6.70 at five loop. The five-loop results with 
imposing ?2"'"(2,2) = 2 and n^(2,2) = 4 are quite close, signaling a weak dependence from 
the exact value of n~^{2, 2) and making the constrained analysis quite safe. 

In order to give a final estimate with a proper error bar, we report in Fig. 1 all our 
results. It appears that all the employed summation methods are slowly converging to the 
same asymptotic value, in particular PBL converges from below, while the other three from 
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TABLE III. Non-defective PBL results for n+(m = 3,3). 



Loop 
Fade 


3 
[1/1] 


[1/2] 


4 


[2/1] 


[1/3] 


5 


[3/1] 


b=0 


6.448 


7.591 




8.152 


8.621 




9.176 


b=l 


6.654 


7.705 




8.134 


8.834 




9.191 


b=2 


6.763 


7.760 




8.123 


8.954 




9.202 


b=3 


6.831 


7.792 




8.116 


9.031 




9.209 



above and quickly. Note in particular that Eq. (8) (method II in figure) is quite fiat up 
to the considered order, suggesting a good convergence. We take as final estimate the very 
conservative one n"'"(2,3) = 6.1(6), which includes all five-loop data. The PBL result is 
only marginally compatible with the final estimate, but this is not surprising since PBL is 
the method that converges slowly and it is expected to underestimate error bars from the 
analysis of n~{2, 3) (note also that the three-loop PBL 3.39 is the almost half of the common 
accepted value). The final value is only marginally compatible with the three-loop estimate 
5.3(2) of Ref. [10]. Furthermore our quoted error is three times bigger than the three-loop 
one. We beheve that our estimate is safe, contrarily to Ref. [10] were probably the valuation 
of the error bar was too optimistic. 

The problem of the value of n"'"(2, 2) may be also exploited resumming the series for 
at e = 2, assuming a smooth dependence on ci up to e = 2. The behavior of the series is 
worst than that for e = 1, in fact Eqs. (6) and (7) are not effective anymore. However we can 
obtain useful informations from PBL method. The approximants [1/1], [1/2], and [1/3] give 
negative values and must be discarded. Only two working approximants remain, namely the 
[2/1] which leads to n+(2,2) = 3.0 and the [3/1] providing n^(2,2) = 3.3, independently 
from b. As a consequence we have a weak indication for n+(2, 2) greater than the conjectured 
value 2. 

B. Nonplanar ordering: m=3 

Let us now consider the case m = 3 which is relevant for spin systems which possess 
nonplanar orderings of the ground state. One possible example of such noncoplanar critical- 
ity was supposed to be the pyrochlore antiferromagnets FcFs [36] and [1] (page 4733). The 
coefficients nf{3) arc (the analytic expressions are in the Appendix C) 

n+(3, 4 - e) = 32.4919 - 33.7184e + 11.1002e^ - 2.1440e^ + 5.2756e^ + O(e^), (9) 
n-(3, 4 - e) = 1.5081 - 0.2816e + 0.5827e2 - 1.4192e2 + 4.0193e^ + O(e^). 

Again for e — 1, PBL technique for n~(3, 3) seems very stable, leading to n~(3, 3) — 1.409(1). 
For n^(3, 3) wc can repeat the same analyses used for m — 2. The PBL results are reported 
in Tab. III. In this case they are much more sensitive to the used approximant and a reliable 
estimate is difficult. Anyway this analysis suggests n+ ~ 9. 
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As for m = 2, the situation is better with considering 
1 



n+(3,4-e) 



0.0308 + 0.0319 e + 0.0226 1 + 0.01460 + 0.0045 e + O(e'), (10) 



which gives rise to n+(3, 3) = 11.7 at three-loop, n"'"(3, 3) = 10.0 at four-loop and n+(3, 3) = 
9.57 at five-loop. 

Following Ref. [10] we impose in the constrained analysis only n+(3,2) = 2, since (at 
least for what we known) this value was never criticized in literature. We obtain 

n+(3, 4 - e) = 2 + (2 - e) (15.2460 - 9.2362 e + 0.9320 - 0.6060 + 2.3348 e^) + O(e^) 
= 2 + (2 - e) a(m = 3, e) , (11) 



and again 

1 



a(3, e) 



= 0.0656 + 0.0397e + 0.0201e^ + 0.0123e^ - 0.0022e^ + O(e^), (12) 



has a better behavior, giving at e = 1, n"'"(3, 3) = 9.975, n"'"(3, 3) = 9.26, n"'"(3, 3) = 9.38 at 
three-, four-, and five-loop respectively. 

The data are much more scattered than those for m = 2. A possible estimate, which 
includes the majority of the previously obtained results, is n"'"(3,3) ~ 9.5(5), that is lower 
than the six-loop fixed-dimension prediction n"'"(3, 3) = 11.1(6) [12]. 

In order to give a complete picture of the critical behavior for all values of n and m, we 
calculate the expansion of n+(m, 4 — e) m = 4, 5, obtaining (applying the same summation 
methods of before) n"''(4, 3) = 12.7(7) and n"'"(5,3) = 15.7(1.0), that are the values reported 
in Table I. These results and their comparison with fixed dimension ones will be commented 
in the Conclusions. 



III. PSEUDO e EXPANSION FOR n+(m,3) 

In this section we re-analyze the six-loop three-dimensional series of Refs. [18,12] with 
the pseudo-e expansion trick, since, as we have mentioned in the Introduction, it provided 
very good results for the marginal spin dimensionality in other models. 

The idea behind this trick is very simple: one has only to multiply the linear terms of the 
two /? functions by a parameter r, find the FP's (i.e. the common zeros of the /3's) as series in 
r and analyze the results as in the e expansion. If one is interested in the critical exponents, 
they are obtained as series in r inserting the FP expansion in the appropriate RG functions. 
With this trick the cumulation of the errors coming from the non-exact knowledge of the FP 
and from the uncertainty in the resummation of the exponents is avoided, obtaining very 
precise results. Note that now, differently from e expansion, only the value at r = 1 makes 
sense. The only drawback of this trick is that it is unable to capture physical features of 
the systems that appear at higher perturbative orders. In fact at one-loop the series are 
equal (apart a normalization) to those of the e expansion. Thus the FP's structure may not 
be different from that of e expansion. All the remarks and the limits of e-expansions apply 
also to pseudo-e one. Explicitly the pseudo-e expansion is able only to evaluate n^(m,3) 
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TABLE IV. Pade table for the series n"'"(2,3). The underlined values come out from Pade 
which possess poles on the positive real axis. The location of the poles nearest to r = 1 is reported 
in brackets. The final estimate is n'^(2,3) = 6.22(12). 





N = 


N = 1 


N = 2 


N = 3 


N = A 


A^ = 5 


M = 





21.798 


6.177 


6.439 


6.289 


6.249 


6.219 


M = 


1 


12.698 


6.435 


6.344 


6.236{3.85} 


6.123{1.31} 




M = 


2 


9.827 


6.290{10.38} 


6.230{3.04} 


6.182{1.74} 






M = 


3 


8.463 


6.247{5.81} 


6.154{1.45} 








M = 


4 


7.695 


6.217(3.87} 










M = 


5 


7.221 













and n/^(m, 3) but not Nc2- For this reason the pseudo-e expansion is unable to recover the 
critical behavior for m = 2 and n = 2,3 coming from fixed dimensional field theory, even if 
it is a re-analysis of the same series. 
For m — 2 we obtain 

n+(2, 3) = 21.7980 - 15.6206 r + 0.2621 - 0.1509 - 0.0392 - 0.0299 + 0{t^) . 

(13) 

At least up to the presented number of loops the series does not behave as asymptotic ones 
with factorial growth of coefficients and alternating signs. So one may apply a simple Pade 
resummation of the series. The results are displayed in Tab. IV. Several approximants have 
poles on the positive real axis (these Pade are the underlined numbers) close to r = 1 and 
thus the estimates of n'^ on their basis should be considered unreliable. Anyway some of these 
defective approximants have poles "far" from r = 1, where the series must be evaluated. 
Thus one may expect the presence of a pole not to influence the approximant at r = 1. 
Indeed all such Pade results are very close since lower order. We choose as final estimate the 
average of those six-loop order Pade without poles at r < 2 (excluding those with N — 0, 
giving unreliable results) , and as error bar we take the maximum deviation from the average 
of four- and five-loop Pade. Within this procedure we have n~^{2, 3) = 6.22(12). In order to 
corroborate this value we resum the series even with the PBL. In this case only the three- 
loop [1,1] and the four-loop [2,1] approximants are nondefective, giving n"'"(2,3) = 6.43 and 
n"'"(2, 3) = 6.35 respectively, in agreement with the Pade estimate at the same perturbative 
order. Anyway a six-loop estimate is not possible within PBL. 

To find n~^{m, 3) for systems with nonplanar spin orderings, we apply the same procedure 
exploited for m = 2 obtaining the following perturbative expressions up to 0{t^) 

n+(3, 3) = 32.4919 - 22.4789 r + 0.4633t^ - 0.2628t^ - 0.1429 - 0.1123T^ (14) 

Again, to extract quantitative information from them, we apply the procedure of before and 
check the results with PBL method. The final estimates are displayed in Tab. I together 
with those coming from five-loop e expansion and other field theoretical methods, also for 
m = 4, 5. 
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IV. CONCLUSIONS 



We compute the Renormalization Group functions of a Landau- Ginzburg- Wilson Hamil- 
tonian with 0(n)xO(m) symmetry up to five-loop in Minimal Subtraction scheme. Such 
higher order computation allowed us to reconstruct with good precision the line n'^{m,d) 
in e expansion up to O(e^). We have also re- analyzed six-loop fixed dimension series with 
pseudo-e expansion. All the results for n~^{m,3) with m = 2,3,4,5 are reported in Table 
I. The results coming from the 1/n expansion are not so near in magnitude to the e and 
pseudo-e ones. Anyway, in this case n+ is determined by perturbative series which are 
known only to 0(l/n^), the correction to the leading term is not small, and since n"*" is 
not expected to be large in magnitude, one is extrapolating the results to value of n which 
could be dangerous. Also the agreement with previous three-loop e-expansion is not good. 
Anyway, as already stated in the text, we believe that the error quoted in previous three 
loop works Refs. [15,10] was underestimated. Differently the agreement is good between all 
the higher order perturbative methods, i.e. direct six-loop fixed dimension, e, and pseudo-e 
expansion. Note that the pseudo-e results are between the five-loop e expansion and direct 
six-loop fixed dimension ones and are the more credible since they are the more stable with 
changing the perturbative order. 

For the physically interesting case of noncoUinear but planar orderings (m = 2), we 
obtain n~^{2,3) = 6.1(6) from five-loop e expansion and n~^{2,3) = 6.22(12) from six-loop 
pseudo-e expansion. Note that the last value seems to exclude n = 6 from the second-order 
phase transition region, in contradiction with Monte Carlo results [17] and EAAM [27] . 
Anyway two remarks on this point are necessary. First, being n"'"(2, 3) — 6 = 0.22(12) small, 
the concept of pseudo-scahng discussed in Ref. [26,4] apphes and the transition is expected 
to be extremely weak first order, i.e. all RG trajectories are attracted toward a small domain 
where the fiow is very slow. In such domain the system spends the majority of its RG time 
(making the correlation length very large) and scaling is partially recovered. Second, the 
estimate of Nc2 — 5.7(3) is marginally compatible with 6, leaving the possibility of a second- 
order phase transition at n = 6 in the scenario of Refs. [18-20]. Anyhow we should expect 
measured critical exponents for n = 6 close to those at n"'"(2,3). Monte Carlo [17] and 
EAAM [27] provide u = 0.700(11), 7 = 1.383(36), and u ~ 0.707, 7 ~ 1.377 respectively, 
that are quite close, but definitely different from the FT perturbative results for n = n'^ 
i'n=n+ — 0.635(4) (see Appendix A). We believe that this apparent disagreement deserves 
further investigations. 
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TABLE V. Critical exponents for n > n"'"(2, 3). Values obtained with other FT methods are 
reported for comparison. The EAAM and MC results for n = 6 are also reported 



1/ 


n = G 




n = 7 


n = 8 


n = 16 


n = 32 


This work 




0.635(4) 


0.71(4) 


0.75(4) 


0.89(4) 


0.94(2) 


6L FD [20] 






0.68(2) 


0.71(1) 


0.863(4) 


0.936(1) 


0(l/n2) [10] 






0.697 


0.743 


0.885 


0.946 


EAAM 


0.707 












MC 


0.700(11) 












7 


n = 6 


T 

n = 


n = 7 


71 = 8 


n = 16 


n = 32 


This work 




1.25(2) 


1.39(6) 


1.45(6) 


1.75(4) 


1.87(4) 


6L FD [20] 






1.31(5) 


1.40(2) 


1.70(1) 


1.860(5) 


0{l/n2) [10] 






1.36 


1.45 


1.75 


1.88 


EAAM 


1.377 












MC 


1.383(36) 












[<^l, W2] 


n = 6 


n = n+ 


n = 7 


n = 8 


n= 16 


n = 32 


This work 




[0,0.86(3)] 


[0.84(3),0.33(10)] 


[0.84(3),0.45(8)] 


[0.86(1), 0.77(2)] 


[0.91(2),0.90(1)] 


6L FD [20] 






[0.83(2),0.23(5)] 


[0.83(2),0.36(4)] 


[0.876(4),0.714(9)] 


[0.933(2),0.868(2)] 


0(l/n) [11] 






[0.768, 0.537] 


[0.797,0.594] 


[0.899,0.797] 


[0.949,0.899] 



APPENDIX A: CRITICAL EXPONENTS 

The critical exponents are computed by expanding in power of e the exponent series at 
the FP. Such a computation gives the exponents only for 

n > (m) = 5m + 2 + 2^6(m + 2)(m - 1), (Al) 

or 

n < no (m) = 5m + 2 - 2^6(m + 2)(m - 1). (A2) 

Indeed, if these bounds are not satisfied the FP's are complex and therefore also the exponent 
series. This is actually only an artifact of the e expansion and the exponents are real 
and well defined for n > n"'"(m, d) (see for details Refs. [7,10]). In order to obtain series 
for the exponents in all the relevant domain we can perform the following trick [10]. For 
< n < Hq we set n = n^{m, 4 — e) + An and re-expand all series in powers of e keeping 
An fixed. In particular, for An = we obtain the critical exponents for n = n+. In the 
most interesting case of m = 2 the exponents along the line n — n'^ are (their analytic forms 
are in the App. C) 

77„=„+ = 0.0208 + 0.0174 + 0.0061 + 0.0032 + O(e^) , (A3) 

and 

(z/„=„+)-i = 2 - 0.5 e + 0.0290 + 0.0767 - 0.0478 + 0.0762 + O(e^) . (A4) 

After a PBL resummation we have z/„=„+ = 0.635(4) and i]n=n+ = 0.049(1). In the same way 
we have also estimated the subleading exponents governing the corrections to the scaling (i.e. 
the eigenvalues of the matrix (B4)). The smallest one is obviously a;"""^ = 0, instead the 
biggest one is (^g""'"' = 0.86(3), after PBL resummation. Now we may calculate critical 
exponents for all n > n+(2,3), even if n < no"(2) ~ 21.8. The resulting critical exponents 
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have a double source of error: one coming from uncertainty in the resummation and one 
from the not precise knowledge of n"'"(2, 3). Note that the latter dominates the given error 
for n — 7, 8, 16, since the resummation of the exponents is very precise, contrarily to that 
of n'^{2,3). To give an example, at n = 7 we have An = 0.9(6), that leads to ^(An — 
0.3) ~ 0.67 and u{An = 1.5) ~ 0.75, and thus the resulting estimate z/(n = 7) = 0.71(4). 
All other exponents are calculated in the same way. For n > n([(2) a direct resummation of 
the critical exponents is possible. In Table V we report the results for n = 7, 8, 16, 32 and 
the corresponding six- loop fixed-dimension and 1/n ones for comparison. There is a quite 
good agreement between all FT methods, a part from the estimates of a; at n = 16, 32, that 
seems to differ significantly from fixed dimension results. The source of such disagreement 
is not completely clear to us. It may be due to an underestimating of the error (especially in 
fixed dimension), or to the bad behavior of the resummation. Anyway it is not so surprising, 
since already for the 0{n) model the estimates of cu are not very good (see e.g. Ref. [3]). 



APPENDIX B: PERTURB ATI VE SERIES 



We calculate the perturbative RG function in the minimal subtraction (MS) renormal- 
ization scheme for the massless theory. We compute the divergent part of the irreducible 
two-point functions of the field (p, of the two-point correlation functions with insertions of 
the quadratic operators 0^, and of the two independent four-point correlation functions. The 
diagrams contributing to this calculation are 162 for the four-point functions and 26 for the 
two-point one. We handle them with a symbolic manipulation program, which generate the 
diagrams and compute the symmetry and group factors of each of them. We use the results 
of Ref. [37], where the primitive divergent parts of all integrals appearing in our computation 
are reported. We determine the renormalization constant associated with the fields (p, 
the renormalization constant Zt of the quadratic operator 0^, and the renormalized quartic 
couplings u,v. The functions (3^, r]fj, and r]t are determined using the relations 

(Bl) 



/3u{u,v) = ^ 



Itn .Vc\ 



(3v{u,v) = ^ 



9 log// 



Vt{u,v) 



uo,vo 



dlogZt 
dlogfj, 



(B2) 



uo,vo 



The zeroes {u*,v*) of the /3 functions provide the FP's of the theory. In the framework 
of the e expansion, they are obtained as perturbative expansions in e and then are inserted 
in the RG functions to determine the e expansion of the critical exponents: 



r] ^r]4,(u*,v*), u ^ (2 - r]4,(u*,v*) +r]t(u*,v*)) ^ 
The stability of each FP is controlled by the 2x2 matrix 

du dv 



(B3) 



d(3y{u,v) df3y{u,v) 



V du 



(B4) 



dv ' 
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A stable FP must have two eigenvalues with positive real parts, while a FP possessing 
eigenvalues with nonvanishing imaginary parts is called focus. 

In reporting the series we make explicit use of the symmetry under the exchange (m, n) 
and write the RG functions as 



= + + + P2 + /Ss" + PI = E u'v'-\Uk-i,i{m, n) + U^-iM m)). (B5) 

(5, = (51 + (51i + (51 + (51 + (51 + (51 (51 = u'v'-\Vu-i,i{m. n) + Vu-iM m)). (B6) 



i=0 
k 



Vcp ^ Vi + V2 + + + rj5, Vk^ Y,'^''"'' \Ek-i,i{m, n) + Ek-i,i{n, m)), (B7) 
1 



1=0 

k 

= i/i + i/2 + i/3 + 1^4 + 1^5, i^fe = Y i^k-i,i{'>^^ ^) + Bk-i,i{n, m)), (B8) 

i=0 



where the first three orders coincide with those of Refs. [7,6,15,10], and the four and five-loop 
coefficients are reported in Tabs. VI, VII, VIII, IX, X.^ 

In order to verify the exactness of our perturbative series we perform several checks on 
them. They reduce to the existing O(e^) ones for the 0(A'^)-symmetric theory [38] in the 
proper limit. For the particular case n = m = 2 they agree, according to the exact mapping 
of Ref. [16], with the four-loop e expansion of the so called tetragonal model Ref. [39]. The 
value of n//(m, 4 — e) coincides with Nc/m, where N^. is the marginal spin dimensionality of 
the cubic model obtained in [40], according to the symmetry argument of Refs. [8,9]. Finally, 
the 1/n expansion of the series coincides with the expansion close to four dimensions of the 
1/n exponents of Refs. [10,11]. 

Since the e expansion is asymptotic, the series must be properly resummed to provide re- 
sults for three-dimensional systems. The first method that is applied to those series that do 
not yet behave as asymptotic is Pade summation, which consists in extending the series ana- 
lytically by the ratio of two polynomials. The approximant [N/M]{t) = Y^iLo ch'T^ / S^lo ^j'^^ 
is determined by imposing the equivalence between its Taylor coefficients with the ones of 
the original series. We also use Pade-Borel-Leroy (PBL) resummation that works as follows. 
Let S{t) = J2 Sk^^ be the quantity we want to resum. The Borel-Leroy transform of S{e), 
is defined by Bij{t) — J2k Skt^ /^{k -1-6-1-1). In order to evaluate S{e— 1) it is necessary to 
perform an analytic continuation of B},{t), that may be achieved using Pade approximants 
[N/M]. Named Ps{b, N, M){t) such analytic continuation, estimates of S{e) are given by 



i 



oo 

dtt''e-'Ps{b,N,M){et) , (B9) 



depending on the considered Pade [N/M] and on the value of the free parameter b. Since the 
integral (B9) should be defined, the approximant must not have poles on the real positive 
axis. The Pade with real positive poles are called defective and must be discarded in the 
average procedure. 



^The complete list of the series is available in electronic format on request. 
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APPENDIX C: ANALYTIC EXPRESSION 



In this appendix we report some useful analytic expressions of the quantities reported in 
the text. We believe that analytical results are very useful for those who want to check our 
results or to compare with other methods where analytic limits are known. 

The analytic expressions of Eqs. (5) are 

n^(2) = 12±4v^, (CI) 

n±(2) = -12Ty%/6, (C2) 

274±91V6 (156±47V6) 
2V y 60 ^ ^ ^ 

± 72771 ± 23369 V6 (156±47\/6) 4 (15004 ± 10621 x/6) 
''^^^>~ 67500 + 7200 + 18000 ^^^^ 

^^^^^^°^^C(5), (C4) 



90 

± _ 61636782 ± 14010773 ^6 (15004 ± 10621 ^6) 4 (464 ± 107 \/6) g 
^"^^ >~ 162000000 ^ 2160000 ^ 68040 ^ 

86358388 ± 29909737 ^6 , . ^ (84368 ± 38339 ^6) 



21600000 ' 720000 
, (2753732 ± 735873 x/6) . , , , 49 (1332 ± 119 x/6) . ,^ , 

+ 432000 ^^^^ + 6400 ^^^^ 

and for Eqs. (9) 

n±(3) = 17±4v^, (C6) 



4(3). -ir.?^. (C7) 



±(o) _ 2267955 ± 461534 V15 42741 ± 10183 V15 , . 

2323200 ^ 10560 ^' ^ ^ 



, 141926963025 ± 31936343543 V15 42741 ± 10183 V15 4 

n?(3) = 1 TT 

^ ^ ^ 101198592000 1267200 



184976182 ± 103127153 V15 _ 4835795 ± 1072837 V15 

^ 306662400 696960 ^ ^ ^ 



17531276380245840 ± 3847065393512279 V15 4835795 ± 1072837 V15 « 
^ 47020913786880000 526901760 



184976182 ± 103127153 V15 4 _ 41214520769938 ± 9458643900547 V15 . . 
^ 36799488000 ^ 7124380876800 ' 



1871364718325 ± 375167813937 V15 , . 49(210903545 ± 35255789 V15) , . 
^ 161917747200 1090355200 ' 



739213610260 ± 238830577907 V15 ,,^,2 

H — C(3) . (CIO) 

539725824000 ^ ^ ^ ^ 

Similarly for (A3), i.e. the exponents for m = 2 and n = n+(2, 4 — e), 
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1 2 5 o /-119 + 24^/6,,^, 3691 + 124 v/6\ , 
48 288 \ 14400 ^ ^ 360000 y 

/6408547 + 408608^6 -119 + 24^ 4 
^ V 1296000000 ^ 1728000 ^ 

-324561 + 51976 \/6,,„, , 5739 - 1508 \/6 5 g 



C(3) + C(5) + O(e^) , (Cll) 



32400000 ' 259200 



and for (A4) 



1 ^ e -l + \/6 2 /37 + 4V6^,,, 397- 37 V6\ 3 / 1104341 - 315351 V6 



18000000 



37 + 4x/6 4 -112429 + 46574 2473 + 484 ... A 4 

^ 120000 ^ ^ 1800000 ^ 36000 ^ 

/19528133- 6267483^6 -112429 + 46574^6 4 2473 + 484^6 g 
^\ 270000000 ^ 216000000 ^ 27216000 ^ 

31404323- 14177332 V6,,„, , 76622 - 74335 a/6 2 



C(3) + ....... C(3)- 



216000000 ' 4500000 



4966979 - 19 1916 1 a/6 49 135 + 7 W6 

21600000 ' 80000 ^ ^ ' ^ ^ ^ ^ 
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TABLE VI. Four-loop coefficients for the RG functions and /S, 



k,i 


Ui^ i{m, n) 


5,0 


(17027 - 51250m + 17042m^ + 154m^ + 38908mn - 26282m^ n - 284 n + 4509 + 222 - 46 n^)/62208 
+7r"'(133 - 436 m + 170 + 331 m n - 226 n + 28 m? n2)/311040 

^C(3)(921 - 2510 m + 711 rn^ - 16 rn'^ + 1763 ni ri - 1096 rri^ n + SO m"' n + 214 r,'-' - 16 rri'' ii'^ + 6 rn'^ /r'')/51Sl 
-C,(5)(1080 - 2820 III + 6()0 in- + 189.') ni n - 970 nr n + loo ;i-)/r)181 


4 1 


(' — '=14n94Q 4- Q64J^9J^ m — 9J^'^n^9 rrj'^ — 1 978 ryt^^ - 7n'^n79 m r) 4- 4^^891 8 m-^ « 4- '^CiQR m-^ n — 7^Q97 m'^ rt'^ 

-3312 + 746m3 n3)/248832 

+7r''(37 + m - 75 - 45 m n + 89 n - 7™^ n2)/77760 

+C(3)(-5511 + 14024m - 3088 + 86 - 9401 mn + 4968 n - lOOm'^n - lOiam^n^ + lb2m^v? - 24m3n3)/5184 
+(^^5^^ — 3615 + 8790m — 1560m2 — 5575m n + 2360m2 n — 400m-^ r7-^V2592 


3,2 


(92919 - 234952m + 49028m^ + Sem''' + 159500mn - 83697m^n - 351m^n + 17059m^n^ + 551m^n^ - 143m^n^)/31104 
+7r*(-143 + 294m - Sm^ - 163mn + 32m2n - 12m2n2)/38880 

+C(3) (6404 - 15260 m + 2472 rr? - 20 m^ + 9784 m n - 4381 m^ ji + 73 m^ n + 999 ir? - 89 'w? + 18 n^) /2592 
+C(5)(1225 — 2750 m + 300 rri^ + 1625 mn — 500 m'^ n + 100 n'^)/324 


2,3 


(-214362 + 477908m - 49184m^ - 296261 mn + 114404 n + 210m^n - 32306m^ - 608m^n^ + 199 n^)/62208 
+7r*(114 - 244m + 16 m^ + 143 m n - 42 m^ n + 13 m^ n2)/19440 

+C(3)(-3854 + 8452 m - 744 m^ - 5096 m n + 1749 rt - 9m3n - 513 + 21 -6m3n3)/1296 
-)-^(5)(_1585 + 3330 m — 160 m'^ — 1845 mn + 360 m^ n — 100 m? n-^)/324 


1,4 


(1 - 2m + mn)((31934 + 3957mn - 14 m^ n^)/15552 - 7r''(170 + 19mn)/38880 
+C(3)(4856 + 552 m n + 3 m^ n2)/2592 + C(5)(480 + 25 m n)/162) 


0,5 


(-196648 - 80456 mn - 6320 + 5m'' n^)/124416 + 7r''(176 + 62mn + 5m^ n^)/38880 

+C(3)(-2332 - 764 rn ii - 63 rii^ ri2)/1296 + C(5)(-930 - 275 ni ii - 10 ni^ ri2)/324 


k,i 




5,0 


(-646929 + 892812m - 105172m^ - 86m'^ - 321700mn + 63034m^n - 94m'^n + 189m^n^ 
4-640m3n2 _ 174m3n3)/248832 

4-7r'*(1409 - 1886 m + 286 + 731 m n - 275 n + 32 n2)/155520 

+C(3)(-18349 + 20756 m- 2402 - 12 - 6203mn + 1720 m^ n + 36m3 n - 166 m^n^ - 36m3 + 6m3 n3)/5184 
4-C(5) (-32250 + 35760m - 1920m2 - 9975mn + 980m2 n - 5m2 rt2)/5184 


4,1 


(275575 - 341500 m + 28538 m^ + 12 m-' + 120765 mn- 18548 m^ n + 22m^n + 89m^ - 240m'' + 87m'' n'')/31104 

+7r'*(-5004 + 5732 m - 648 m^ - 2147 m n + 784 m^ n - 112 m^ 11^)1X5552^ 

+C(3) (32300 - 32284 m + 2730 m^ 4- 6 m^ 4- 9321 m n - 2380 m^ n - 36 m^ n 4- 312 tr? 4- 54 m^ r? - 12 m^ n3)/2592 
4-C(5) (54940 - 55940 m 4- 2160 m^ 4- 15225 mn - 1160 m^ n)/2592 


3,2 


(-797194 + 778564 m - 34544 m^ - 278829 mn + 29972 rn^ n - 22 m'' n - 728 m^ 4-480m''n^ - 261 m'' n'')/62208 
4-vr*(1806 - 1660m4- 96m2 + 607mn - 194 m^ n 4- 38m2 ri2)/38880 

4-C(3)(-23278 + 18624 m - 836 m^ - 5361 mn 4- 1149 m^ n 4- 9 m^ n - 222 m^ - 27m3 4- 9m3 n3)/1296 

4-C(5)(-9545 + 8010 m - 160 m^ - 2190 m n 4- 120 m^ n)/324 


2,3 


(294788 - 179096 m + 75554 m Ti - 4808 m^ n + 269 m^ -SOm" 4- 87m^ n^)/31104 
4-vr*(-340 4- 190m - 86 mn 4- 17 m^ n - ern'^ n2)/9720 

4-C(3) (17368 - 8416 m + 2970 m n - 402 m^ n 4- 144 m^ n^ 4- 9m3 n^ - 6m3 n3)/i296 
4-C(5)(1735 - 880 m + 295 m n - 10m2n)/81 


1,4 


(-157432 - 27832 mn - 56 m^ n^ - 29 m'' n'')/41472 4- 7r"'(90 + 17 rnn + m^ n^)/6480 
4-C(3)(-4536 - 560mn- 24m2 n^ 4- m^ n=')/864 - C(5)(910 4- 105mn)/108 
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TABLE VII. Five-loop coefficients for the RG function 



k,i 


Uk,i(m, n) 


6,0 


(7720646 - 26391892m + 11963785m^ - 1016454m^ + 3269m* + 22679912m n - 20725184 m^n + 1764720 n 

— 1 4Do (71 77. -|- 4 / o'iiOo 772 77 — oiDUy^ 777 77 -|- yZiO TTL TL -\- O ( Li)Z Til Ti — OOOo 777 77 -|- l 1 1 Til 77 if iiy^loyoO 

+7r'l (6643 - 19108 m + 3786 m? + 2284 - 248 m* + 10402 m n + 3006 rn^n- 5222 n + 520 m" n - 5395 'w? v? 
+4462 m'^n^ - 464 m'' - 882 m^n^ + 240 m"' - 24m* n'')/i4929920 

+7r''(2705 - 9140m + 3520 m^ + 210 m"^ + 6950 m n - 4155 m'^ n ~ m5m'^n + 40m'^ + 555m'^ - 80m^ n3)/15676416 
+C(3)(180862 - 544185 m + 204185 m^ - 21939 m^ + 215 m"' + 417829 mn - 332055 m^ n + 40996 m^ n - 414m'* n 
+76436 m^ n^ - 25305 m'^ n^ + 303 m" n^ + 3187 m^ - 126 m* n^ + 11 m" n"')/248832 

+C(3)2(-37 + 1048m - 1292m2 + SlSm^ - 1516m n + 2571m2 n - 587m3 n - 822m2 n^ +365m3 v? - 48 m^ n3)/41472 
+C(5) (269531 - 852440 m + 353290m2 - 39912m-^ + 701125 mn - 627204 n + 77394 n 
+160298 m^ n^ - 46682 m'^ n^ + 4600 m'^ n3)/248832 

+C(7)(79331 - 228438 m + 69776 m2 + 167139 mn - 105840 m^ n + 18032 m^ n2)/55296 


5,1 


(-23038930 + 73796176 m - 29373553 m^ + 1658632 m^ - 3395 m* - 60585772 m n + 50569562 m^ n - 3204600 m-^ n 
+ 10406 m*n - 11479472 m^ n^ + 1780134 m^ n^ - 17199 m^n^ - 124124 m^ n'^ + 14082 m" n'^ - 1947 m" n4)/5971968 

+71"' (23753 - 63434 m + 21402 m^ - 5710 m^ + 236 m* + 55489 m n - 60768 m^ n + 13876 m^ n - 652 m* n + 25316 n^ 
"1 9nfifi tfi3 ,.,2 _i_ Qfin ■)-i-i4 ,,.,2 _i_ oono tfi3 .^,3 f;n/i i-vi4 ^3 i ro '>t-i4 ^4\ /'z/ii^/iQRn 

— iZUOO Afi /t -\- III II -\- Z.Z\)Z ilL il — OU'4: //i ll DU A/t M / '4:D'4:yDU 

+7r*' (16910 - 23540 m - 5480 in^ - 4800 m'^ + 13575 m n - 15650 m^ n + 12040 m'^ n + 15625 m^ n^ 
-10120 m-^ n^ + 1440 m'^ n-^)/47029248 

+C(3)(- 1176575 + 3390386 m - 1109336 m^ + 72418 m» - 318 m* - 2523565 m n + 1806256 n - 150316 m^ n 
+804 fre" n - 402062 m^ n^ + 108034 m^ n^ - 830 m" n^ - 15286 m'' n'' + 436 m"* n'' - 46 m"* n*)/248832 
+C(3)2 (-28978 + 57820m + 4744 m^ -4608m3 - 24477mn - 17978 m^ n + 9112 m^ n 
+9733 m2 n2 - 6232 + 864 m^ n3)/124416 

+C(5)(-406105 + 1201928 m- 432714 +42996 - 947945 m n + 783954 n - 89992 n - 206134 
+61028 m^n^ - 7016m3n3)/62208 

+C(7)(-44786 + 124852 m - 35280 - 89621 m n + 54390m2n - 9555 n2)/4608 


4,2 


(115482976 - 343572912 m + 116333720 m^ - 3730368 + 3608m* + 268202738 mn - 201574416 n + 8761141 m'^n 
-19289 m* n + 45741288 n^ - 6288827 + 46947 m* n^ + 652565 ■n? - 47076 m* + 7905 m* n*) /11943936 
+7r4(-267336 + 708928 m- 194576 + 20544 - 224 m* - 551881 m n + 449912 m^ n - 56150 m^ n + 1072 m* n 

— iOv54zy 777 77 -|- OOOZZ 777 77 — ZUUU 777 77 — yM4l:?77 77 -f- iOOZ 777 77 — Z4U 777 77 ^/l^yzyazU 

+7rS (-174840 + 412360 m - 71640 + 8960 - 281095 m n + 174430 m^ n - 24600 m^ n - 62815 m^ n^ 
+22840 m^ - 3600 m^ n3)/47029248 

+C(3)(6054444 - 16659432 m + 4720608 - 170240 m^ + 176 m* + 12109993 mn- 7988808 m^ n + 429072 n 
-818 m* n + 1824513 m2 n^ - 382176 n^ + 1350 m^n^ + 62156 m^ n^ - 968 m^n^ + 130 m* n*)/497664 
+C(3)2 (98640 - 215480 m + 11160 m^ + 7040 m^ + 111485 mn + 8998 m^ n - 16488 m^ n - 16963 m^ n^ 
+13768 m^ n^ - 2160 m'^ n-^)/124416 

+C(5)(1143847 - 3112678 m + 887044 m^ - 62060 m» + 2306576 mn - 1653546 m^ n + 153072 m^ n + 446461 m^ n^ 
-126420 m^ n^ + 17704 m'' n'^)/62208 

+C(7)(125189 - 331338 m + 80360 m^ + 229957 mn - 128576 m^ n + 24108 m^ n2)/4608 


3,3 


(-41026088 + 109185000 m - 27504368 m^ + 371544 m^ - 79647772 m n + 51460526 m^ n - 1351466 m^ n + 1484 m* n 
-12678205 m^n^ + 1407142 m=^ n^ - 6890 m*n2 - 218338 m'^ n^ + 9456 m* n^ - 2025 m* n*)/2985984 

H-TT I^OZy 14 — IDiOoU 771 H- v5 ( v5^o ?77 — 10(D?77 -|- iZUoD ( TflTl — oOvSZD TTt 77 -|- OZUO 777 77 — OO 777 77-i-ZMol?77 77 

-7596 m^n^ + 132 m* n^ + 1560 n^ - 156 m*n3 + 30 m* n*)/1866240 

+7r6 (23335 - 55490m + 9140m2 - 320 m^ + 37255 mn - 20300 m^n + 1280 m^ n + 6360 m^ - 1560 m^ n^ 

+300m-^n-^)/2939328 

+C(3)(-1055134 + 2705500 m - 604488 m^ + 9256 m^ - 1901069 mn + 1132320 m^ n - 35686 m^ n + 4 m* n 
-285779 m^ n^ + 43744 m^ n^ - 18 m* n^ - 8669 m^ n^ + 24 m* n'* - 5 m* n*) /62208 

-C(3)2(10409 - 22174m + 1132m2 + 224m3 + 11441mn + 92m2n - SOOm'^n - 1080m2n2 + 'imm^ri^ - ISOm^n^) /777& 
+^(5)(-454646 + 1130856 m - 229244 m^ +7680m^ - 781415 m n + 460582 m^ n - 28608 m^n - 132069 n^ 
+32800 ni'' n'^ - 5936 m^ n»)/15552 

+C(7)(-36946 + 89572 m - 15680 m^ - 58653 m n + 27734 m^ n - 6027 m^ n2)/864 


2,4 


(140915584 - 325697728 m + 43866560 m^ + 217044656 mn - 109698320 n + 1306736 m" n + 34192982 m^ 
— zoouuuu m n +o/yom n +Doo/(om n — i4zoDm n +4z4om n )/iiy4oyoD 

+7r* (-65666 + 153224 m - 21892 m^ - 106391 m n + 61604 n - 2046 m^ n - 22040 m^ n^ + 4392 m^ - 24 m* n^ 
-1200 m^ + 54 m* n^ - 15 m* n*)/1866240 

+7r'5(-8245 + 18410 m - 1920 m^ - 11770 mn + 5260 m^ n - 130 n - 1810 m^ + 280 m^ - 75 m^ n3)/979776 

-rs / I^ODZC5410 — CSolOoOZ 777 -j- lUOyOZU 777 -p 04 / 1UZ4 777 77 — ZDD040D 777 77 -p O / i DU 777 77 -p C54oOo4 777 77 — OOUOD 777 77 

-76 m* + 22575 m^ + 146 m* - 35 m* n*) /248832 

+C(3)2(3409 -6898 m + 80 m^ +3370 mn + 236 m^n- 78 m-^n - 242 n^ + 168 m^ - 45 m^ n3)/2592 
+C(5)(283674 - 637364m + 70016m2 +407131mn- 182346m2n + 5448 m3n + 61874m2 - 11418m3 
+2985 m^ n3)/10368 

+C(7)(45815 - 100450m + 8820 m^ +60760mn- 21070m2n + 6125 m^ n2)/1152 


1,5 


(1 - 2m + mn)((-34395968 - 6592896 mn - 119350 m^n^ - 483 m^ n^)/5971968 

+7r*(40708 + 9356mn + 513m2 n^ + 3m3 n3)/1866240 + 7r'5(4810 + 815 m n + 30 m^ n2)/97977g 
+C(3)(-934072 - 172480 mn - 4080 m^n^ + 17 m'^ n3)/124416 + C(3)^ (-1574 + 95 mn + 18 m^ n2)/2592 

+C(5){-37618 - 6069 mn - 301 m^ n2)/2592 + C(7)(-5782 - 637mn)/288) 


0,6 


(13177344 + 6646336 m n + 808496 m^ n^ + 12578 nr' n'' + 13 m* n*)/3981312 

+7r*(-21120 - 9532mn- 1388 n^ - 63 m^ n3)/1244160 + 7r«(-7440 - 3130 mn - 355 m^ n^ - 10 m^ n3)/1959552 
+C(3)(1314336 + 552280 m n + 67640 m^ n^ + 1248 m^ n» - 9 m* n*)/248832 

+C (3) ^(3264 + 446 mn - 59m2 n'^ -6m^ n3)/5184 + C(5)(165084 + 66986 mn + 7466 m^ n^ +305m3 n3)/15552 
+C(7)(25774 + 9261 mn + 686 m^ n2)/i728 
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TABLE VIII. Five-loop coefficients for the RG function 



k, i 




6,0 


(-72238776 + 132078652m - 26680036m^ + 273681m-^ - 345m"' - 60428904 mn + 22666801m^ n + 204824m^ n 
-533m'* n - 1463603 rn^ - 234959 ^ 4375 „4 „2 _ ^^^2 nv^ - 4626 m"* + 732 m'^ n*)/11943936 
+77"* (283706 - 451827 m + 105356 ,n'^ - 5907 - 24 m"' + 200762 mn - 107414 m^ ji + 6718 n + 96 n 
+17168 m^ - 2650 m^ - 144 + 84 + 96 - 12 n*)/7464960 
+7r6 (435460 - 674770 m + 130820 - 5040 + 288075 m n - 128455 n + 4790 n 
+16985 m^ ri2 _ 1175 n2)/47029248 

+C(3)(-5923864 + 8822466 m - 1755312 m^ + 19168 m'^ +334m4 - 3429421 mn + 1453852 m^ n - 4514 m^n 
-734 m^n- 170929 r? - 6198 m^ 752 + 2384 m'^ n'^ - 388 m"' n'^ + 40 m"* n'*)/497664 

+C(3)2 (-317108 + 284654m - 6772m2 - 2640m-^ - 40197mn - 14179m2 n + 2426m-^ n + 3431m2n2 - 641m3 n2)/124416 
+C(5) (-6417032 + 9071536 m - 1459068 m^ + 81384 m"^ - 3393769 m n + 1212152 n - 67088 n 

1 "^Qfif^fi ^ 2 ^1 2 _|_ 1 Rovn im'^ t> fins ^ ssQO 
— ±ooDDD ni II -\- iDUiU III ii — nvjo III 11 jj Z'ioooz 

+C(7) (-2229108 + 2985472 m - 375144 m^ - 1019788 m n + 264208 n - 18277 n2)/55296 


5,1 


(77055030 - 125871804m + 21294903m^ - 149415m^ + 96mi + 55861516m n - 19092239 m^ n - 180275 m^n-V 262 n 
+1334594 m^ + 213766 m'^ - 3282 m"' rp- - 4687 m^ n'^ + 4626 m* ra^ - 915 m* n4)/2985984 
+7r4(-620803 + 882404 m - 171586 n? + 6676 m» + 12 m* - 381314 m n + 192706 m^ n - 10104 n - 96 n 
-33043 n? + 5296 m^ + 216 m"* - 324 m^ - 192 m"' + 30 m* n"')/3732480 
+7r'^ (-936130 + 1316400 m - 210360 + 5600 - 548055 m n + 230560 n - 6720 n 
-33685 + 2120 n2)/23514624 

+C(3)(6309091 - 8414600 m + 1418328 - 10638 - 88 m* + 3170292 m n - 1282310 n + 2934 n 

+364 n + 168641 + 5428 - 564 m* - 2348 + 388 - 50 n4)/124416 

+C(3)2(619622 - 535776m + 13176m2 + 2912m3 + 65193mn + SOTeOm^n - 3360m3 n - 7237m2 + 1208m3n2)/62208 

+C(5) (3418732 - 4467538 m + 584784 - 22800 + 1633807 mn - 535301 m^ n + 22976 m^ n 

-66871 rn^ - 7352 rn^ + 321 m'^ ri*)/31104 

1 )[l Ut)Ut)0 — lAJoiiU IfL -\- L()L:)Z,c> III -\- t>ioO'J. ( III IL — 1 :A)uo III fi-\-:.il-)'.}ni fl ) / iOUo 


4,2 


(-139217840 + 194720568 m - 24121596 + 91684 m'^ - 83985744 mn + 24543584 n + 200423 m'^ n - 131 m* n 
-2334581 _ 291418 + 3282 + 23336 - 6939 + 1830 n*)/2985984 
+7r"'(2301104 - 2794568 m + 402488 m^ - 7856 + 1178888 mn- 545040 n + 20644 m^ n + 96 n + 104065 n? 
-15916 m3 - 432 m" + 1632 m^ + 576 m"' - 120 n*)/7464960 

+7r6 (3382640 - 4147520 m + 486480 m^ - 6400 m^ + 1684280 ran- 650500 m^ n + 13040 m^ n 
+109785m2 - 5960 n2)/47029248 

+C(3) (-22967160 + 26065200 m - 3253224 m^ + 13408 - 9525746 m n + 3532984 m^ n - 5288 m^ n - 364 m"' n 
-553459 m^ n^ - 15716 m'^ + 1128 m^n^ + 7540 m^m^ - 1164 m^'n'^ + 200 m"' n"')/248832 
+C(3)2 (-2079760 + 1637632 m - 28080 m^ - 3328 m» - 169288 mn - 106228 m^ n + 6416 m^ n 
+24069 m^ n^ -3512m3 n2)/124416 

^^0 j — \.LL^^±^j^\} i^ioiUoD AM — - loOOoOD III ~r ZD'4iOU /fi — OUooZ ( O //t /I ~r i'loU'4iyD III II — ^OO^^J TTl Th 

-219687 m^ n^ + 21140 m-^ n^ - 1128 m'^ n'^)/62208 

+C(7)(-1342404 + 1494108 m - 115640 - 484267 nt n + 98196 m^ n - 9408 m^ n2)/4608 


3,3 


(67986200 - 75382408 m + 5048536 nt^ + 32379740 m n - 7295478 n - 34186 n + 1144333 + 88907 
-547 m* n^ - 16178 m^ n^ + 2313 m"* n^ - 915 m"* n"*) /1492992 

+7r"' (-581572 + 554120 m - 43136 m^ - 236664 m n + 92168 n - 1808 n - 21599 n^ + 2670 m? + 36 n^ 

-462 m=^ - 96 n^ + 30 m* n'*)/1866240 

+7r'*(-207670 + 203340 m - 12800 m^ - 83285 m n + 27270 m^ n - 280 m^n - 5880 nt^n^ + 240 m^ n2)/2939328 
+C(3) (5715252 - 5053816 m + 339728 m^ + 1855772 mn - 572300 m^ n + 428 m^ n + 119493 m^ n^ +2696m3 
— y-yi /(i 11 — loyi iii ti -t- ly^ iii n — ou rii ii j / ozzvo 

+C(3)2 (120242 - 76548 m + 448 m^ + 6379 m n + 5334 m'^ n - 136 m^ n - 1296 m^ n^ + 144 m^ n2)/7776 
+C(5)(2982544 - 2773064 m + 145016 + 1014902 m n - 244024 n + 3780 n + 47575 _ 3432 ^3 ^2 
+254 n3)/15552 

+C(7) (240933 - 216090 m + 8820 + 70560 m n - 11270 n + 1568 n^) /864 


2,4 


(-152904896 + 107048320 m - 53120672 mn + 7504128 m'' n - 2268736 m'' n^ - 82186 m"" n^ + 36434 m"" n^ 
-2313 n^ + 1830 n*)/5971968 

+7r*(340984 - 196700m + 106176mn - 26092m2n + 9902m2 - 675m3n2 + 243 m^n^ + 24m*n^ - 15m4n*)/1866240 
+7r'^ (39720 - 23650 m + 12170 mn- 2535 n + 910 - 20 n2)/979776 

+C(3)(-6552944 + 3545400 m - 1587312 mn + 309792 n - 110060 _ 1472 r„3 „2 _,_ ;^4gg ^3 „3 
-97 m* n^ + 50 n*)/124416 

+C(3)2 (-21632 + 8110 m - 414 m n - 559 n + 206 - 12 n2)/2592 

+C(5)(-282530 + 159690 m- 73301 mn+ 11355 m2 n - 3739 m^ + 148 n^ - 20 n3)/2592 

+C(7)(-44835 + 24500m - 9800 mn + 980 m^ n - 245 in'^ n'^)/288 


1,5 


(24902080 + 6117632 mn + 222434 m^n^ - 3557m-' n'^ - 183 m"* n"')/2985984 

+7r "'(-108664 - 26368 mn - 2110 m^n^ - 48 m'^n'^ + 3 m"' n"')/ 1866240 + 7rS(-6350 - 1475mn- 95m2 n2y489888 
^C(3)(1037132 + 188904 ;/( ri + 11110 in- rr - 151 ni'' ri'' - 5 r|4')/62208 + C(3)^(3322 - 23 rn ri - 23 rn^ ri-)/1296 
-((5) (81)8 18 - 172:->0 111 11 + 795 m- 11- + I ni' n:'')/2r,n2 + (,(7) (11210 + 2205 ni 11 + 19 rn - ;)-)/288 
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TABLE IX. Four and Five-loop coefficients for tlie RG function ry. 



k, I 




4,0 
3,1 
2,2 
1,3 
0,4 


— 5(— 397 + 1100m — 292m-^ — 14m'^ — 760mn + 394m-^n + 2Qrrv^n — 43m''n'' — Idm'^n'' + 2m'^n'^)/165888 
5(-147 + 376m - SOm^ - 2m^ - 253mn + 12Aw?n + Gm^n - l^m^v? - Gm^n^ + 'm?v?)/2Q7m 
-5(-118 + 268m - 32m2 - 171mn + 72m2n + 2m3n - ISm'^n'^ - im^ri^ + m3n3)/13824 
5(1 - 2m + mn)(-100 - 18mn + m2n2)/20736 
-5(2 + mn)(-100 - 18mn + m'^n^)/829U 


5,0 


(174432 - 550224 m + 207810 m^ - 6801 m^ + 351 m"* + 436969 mn- 332893 m^ n + 9998 n 
-819 m-* n + 63604 m? - 2905 m^ + ygo ^4 „2 ^ 49 ^3 „3 _ 390 ^4 „3 ^ 39 ^4 „4) /5971968 
-t-TT* (314 - 908 m + 280 m2 + 659 m n - 410 m^ n + 65 m^ n2 ) /2488320 

+C(3)(-1490 + 5528m- 2320 m^ - 192 m^ - 36 - 4893 mn + 3738 m^ n + 436 n + 84m* n - 511 
-316 - 80 m* + 16 m^ + 40 m"' ft'' - 4 n4)/497664 


4,1 


(-635344 + 1905920m - 649324m^ + 14482m^ - 390m* - 1482454mn + 1083994m^n - 26566m^n + 1560m*n 
-222127m2n2 + 11924™^™^ - 2340m,^n'^ - 700m3n3 + ISeOm^n^ - 195m4n'')/5971968 
+71* (- 163 + 460 m - 134 m2 - 331 m 71 + 202 n - 34 rp- ) /373248 

+C(3)(1289 - 4590 m + 1894 + 108 + 10 m" +404:'lmn - 3154 m^n - 304m3n - 40 m'^ n + 457 
+286 + 60 rn* - 25 m'^ n'^ - 40 m"' + 5 m* n") /124416 


3,2 


(493416 - 1352120m + 369376m^ - 4088m^ + 1014866mn - 689396m^n + 12174m^n - 390m'*n + 163973m^n^ 
-9096m3n2 + inOrn'^n^ + lOgOm^^n^ - 1170m4n^ + 195m4ra'*)/2985984 
+7r4(61 - 160 m + 38 + 112 mn - 64 n + 13 n2)/93312 

+C(3)(-2032 + 6300 m - 2172 m^ - 64m3 - 5357 m n + 4102 m^ n + 292m-^ra + 20m4 n - 736 m^ - 398 m'^ 
-mm*n^ + 55m3 + 60 m^n^ - 10 m* n4)/124416 


2,3 


(-417648 + 995568m - 160272m^ - 704492m?i + 417416m^?i - 4000m^n - 131402m^n^ + 6050m^n^ - 390m*n^ 
-1415m''^n'^ + 780m4n-^ - 195m'ln4)/2985984 

+7r4(-50 + 116 m - lem^ - 77m n + 38 ra - llm^ n2)/93312 

+C(3) (1856 - 4616 m + 904 m^ + 3594 mn- 2472 m^ n - 100 m^n + 649 m^ + 250 m^ + 20 m* v? 

-55 m^ -40 m'* n'^ + 10 rt4)/124416 


1,4 


(1 - 2m + mn)(5(77056 + 22752irtn + 296m^n^ + 39m^n^) + 647r'*(22 + 5mn) 
-240C(3)(184 + 64mn - Qni^-n? + m^n^)) /bSilimS 


0,5 


-(2 + mn)(5(77056 + 22752mn + 296m^n^ + 39m^n-^) + 647r*(22 + 5mn) 
-240C(3)(184 + 64mn - 6m^n2 + m^n^))/29859840 
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TABLE X. Four and five-loop coefficients for the RG function Xjv. 



Bk.i{m,n) 



(112697 - 323980m + 97772ni^ + 814m''' + 237400mn - 150394m^n - 426m^n + 26503m^n^ - 384m^n^ 
-2m3n^)/497664 -|-7r4(163 -460m+ 134m2 + 331 mn - 202 n + 34 n2)/311040 

+C(3)(273 -632m + 136m2 - SOm^ + 373mn - IQGm^ n + 82m3 n + 52 - 44m3 + Gm^ n3)/10368 

(-42647 + 112736m - 27340m2 - 102m3 - 80233m n + 47674 m^ n + 56 m^ n - 10189 m^ + 44 ^3 „^ ^ ^3 „3)/62208 
+7r4(-61 + 160m -38 m2 - 112 mn + 64m2 n - 13 m^ n2)/38880 

+C(3)(-214 + 492m-78m2 + Um^ - 284 mn + 113 m^ n - 37 m^ n - 35 m^ n^ +35 m^n^ - 6m3n3)/2592 
(35118 - 81708m + 11472 m^ + 55091 mn - 28472 m^ n - 2m3 n + 8498 m^ n^ + 4m3 n^ - n^)/iU72 
+7r"'(50 - 116m + 16m2 + 77mn- 38m2n + llm^ n2)/25920 

+C(3)(82 - 184m + 20m2 + 106mn - 34m2 n + 6m3 n + ISm^ n^ - 12 m^ n^ + 3m3 n3)/864 

-(1 - 2m + mn)(155300 + 37890mn - bm^-n? + 7r4(352 + 80mn) + C(3)(16320 + 2400m« + 720m2n2))/311040 

(2 + mn)(155300 + 37890mn - 5m^n^ + 7r"'(352 + 80mn) + C(3)(16320 + 2400mn + 720m^n^))/1244160 

(6716354 - 21727736 m + 8800405 m^ - 507480 m'^ + 2103 m" + 17902846m n - 15072442m'' n + 995884m'' n - 1398 m* n 

+3424054 r? - 575534 m^ r? - 537 m"' + 43670 m^ - 210 m* + 21 n*)/11943936 

+7r'* (24815 - 76670 m + 28406 - 1226 m'^ - 140 m"* + 60313 mn - 46884 m^ n + 2660 m^ n + 268 m* n + 10142 m^ n^ 
-1582 m'^ Ti^ - 224 m"' + 14 m'^n'^ + 120 m^n^ - 12 m* n*)/14929920 

+7r*' (34410 - 105900 m + 39000 rr? ~ 1920 nfi + 82855 m n - 63370 m^ n + 3560 m^ n + 13165 m^ n^ 
-1960 m^ n^ + 160 m^ n^)/94058496 

+C(3)(156293 - 443202 m + 147804 m^ - 17482 m'^ + 294 + 320673 m n - 231192 n + 33756 m^ n - 708 m'ln 
+52756 rrp- - 22810 r? + 686 + 3438 m^ n^ - 340 m'' + 34 m"* n^) /497664 

+C^(3)(-23838 + 70212 m - 22920 m.^ +384™^ - 52157 m n + 34814 m^ n - 712 m^ n 
-6143 m^ rp- + 392 - 32 n=^)/248832 

+C(5)(96 - 671m + 556 - 77 m-^ + 698 mn - 701 n - lA.m'- n - Vim^ + 171 m'^ ra^ - 35 m^ n^)/31104 



<(5)(96 - 671m + 556 m^ - 77 m'^ + 698 mn - 701 m"^ n - 24m-'' n- 13m'' + 171 m'^ - 35m^ n^)/31104 

-24959332 + 76197848 m - 27373258 m^ + 1096198 m^ - 2124 m* - 61187662 mn + 48929662 m^ n - 2753983 m^ n 

, 1797 m-* n - 11819002 m^ + 2081903 m'' - 303 m"' n^ - 212479 m^ + 840 m"' n'* - 105 m"' i&) /11943936 

+7r-* (-89558 + 264608 m - 88316 m^ + 2672 ni'^ + 152 m" - 204363mn + 152024 m^ n - 7410 m^n - 496 m^n 

-35089 m^ n^ + 5766 m'^ n^ + 704 m** n^ - 274 m'^ n^ - 480 m** n'^ + 60 m"* n"*) /14929920 

+77*^ (-31085 + 91700 m - 30610 m^ + 1080 m^ - 70435 mn + 51610 m^n - 2440 m'^ n 

-11380 m^ n^ + 1760 m^ n^ - 200 m^ n3V23514624 

+C(3)(-571026 + 1564720 m - 462600 m^ + 40260 m^ - 328 m* - 1112181 m n + 750388 m^ n - 92100 m^ n + 1354 m" n 
-182311 m^ n^ + 78880 m'^ n^ - 2046 m"* n^ - 14200 m^ + 1360 m" n'^ - 170 n'')/497664 
+(^ (3) (20887 - 59740 m + 18182 m^ - 216 m^ + 43877 mn- 28502 m^ n + 488 m^ n + 5336 m^ n^ 
-352 m3 n^ + 40 m^ n3)/62208 

+C(5)(-736 + 2826 m - 1524m^ + 170m3 - 2516m n + 2084 m^ n + 122 m^ n - 17m^ n^ - 526 m^ n^ + 117m^ n3)/31104 
' I o^ni Qfi/l _ c;/i«QQn«n ™ _l 1 Kfinnon/l ™^ _ QOnR'70 ™ J _l /1 00*71 _ Qiic?!!/!™^*, j_ nQOooift^-i*, _oin™4*. 



(19701864 - 54683960 m + 15600904 m'' - 320672 + 42271184 mn- 31187114 m-' n + 1328916 n - 210 n 
+8513717 m^ n^ - 1441854 + 630 m* n^ + 217120 m^ - 630 m* n^ + 105 m* n*)/59'''1968 

+-7r*(17255 - 46902 m + 12576 m^ - 184 m^ + 35060 mn- 24131 m^ n + 883 m^ n + 30 m* n + 6323 m^ n^ 
-1001 m^n^ - 90 m^n^ + 106 m^ n^ + 90 n^ - 15m4 n*)/1866240 
+-7r« (11995 - 32600 m + 8770 m^ - 160 m^ + 24220 mn- 16420 m^n + 580 n 
+4135 m^ n^ - 620 m^ n^ + 100 m^ n^) /5878656 

+C(3)(218618 - 564220 m + 133328 m^ - 6344 m^ + 390068 m n - 238048 m^ n + 22302 m^ n - 170 m" n + 65749 m^ n^ 
-27288 m^ n^ + 510 m* n^ + 5920 m^ n^ - 510 m* n^ + 85m* n*)/124416 

+C^(3)(-7889 + 20920m - 5174m2 + 32m3 - 14924mn + 9044m2n - 116m3n - 1997m2n2 + 124m^n'^ - 20m3n3)/15552 
+C(5)(1364 -4444m + 1796m^ - 80m3 + 3569mn - 2486m^ n - 208m3 n + 38m^ n^ + 614m3 n^ - 163m3 n3)/31104 



(-16860112 + 40609072 m - 6888848 m^ - 29640128 m n + 19136204 m'-^ n - 465020 m^ n - 6579128 m^ n^ 
+911110 m^ n^ - 210m* n^ - 223255 m'^n'^ + 420 m*n3 - 105 m* n*)/5971968 
+7r*(-14534 + 34544 m - 5476 m^ - 24459 m n + 14676 m'^ n - 302 m^ n - 4923 m^ n^ 
+616 m^ n^ + 30 m* n^ - 127 m'^ n'^ - 60 m* n'^ + 15 m* n*) /1866240 

+7r6(-10110 + 24060 m - 3840 m^ - 16915 m n + 9970 m^ n - 200 m^ n - 3265 m^ n^ + 400 m^ n^ - 100 m^ n3)/5878656 
+C(3)(-178724 + 417584 m - 60136 m^ - 277386 m n + 145188 m^ n - 8000 m^ n - 51406 m^ n^ + 17930 m^ n^ 
-170 m* n2 - 5135 m3 n^ + 340 m*n3 - 85 m* n*)/124416 

+C2(3)(6522 - 15252 m + 2208 m^ + 10313 mn - 5414 m^ n + 40 m^ n + 1643 m^ n^ - 80 m^ n^ + 20 m^ n3)/15552 
+C(5)(-1320 + 3480 m - 840 m^ - 2486 mn + 1376 m^ n + 116m'^ n - 89 m^ n^ - 358 m^ n^ + 121 m^ n^)/31104 
(1 - 2m + mn)(315(3166528 + 1077120 m n + 45254 m^ n^ + 21 m'^n^) + 10087r*(2668 + 816mn + 29m^n^ -3m^n^) 
+32007r«(186 + 55 m n + 2 m^ n^) + 15120C(3)(31848 + 8208 m n + 940 m^ n^ + 17 m'^n'^) 
-241920C^(3)(582 + 145 mn + 2 m^ n^) + 241920C(5)(72 + 14 mn - 5 m^ n^))/752467968 



(2 + mn) (315(3166528 + 1077120 mn + 45254 m^ n^ + 21 m'' n'') + IOOStt" (2668 + 816 m n + 29 m^ n^ - 3m^ n^) 
+32007r'5(186 + 55mn + 2m2 n^) + 15120C(3)(31818 + 8208 m n + 940 m^ n^ + 17 m'St'^) 
-241920C2(3)(582 + 145 mn + 2 n^) + 241920C(5)(72 + 14mn - 5m2 n2))/3762339840 
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